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Abstract 

Previously, in an effective field theory formulated by us, we have carried out parameter- 
free calculations of a large number of low-energy two-nucleon properties. An experiment 
y . at the Institut Laue-Langevin is currently measuring spin-dependent effects in the po- 

(~| , larized np capture process n + p ^ d + 7 at threshold. Noting that spin-dependent 

observables for this reaction are sensitive to terms of chiral orders higher than hitherto 
studied, we extend our effective theory approach to this process and make parameter-free 
predictions on the spin-dependent observables. 
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1 Introduction 

An experiment by Miiller et al |]l| at the Institut Laue-Langevin (ILL) is currently 
measuring spin-dependent effects in the thermal neutron capture process 

n+p ^ d + '^ , (1) 

wherein, as indicated, both the incident neutron and the target proton are polarized. At 
threshold, the initial nuclear state is in either the ^So{T = 1) or "^SiiT = 0) channel, where 
T is the isospin. The process (ffl) therefore receives contributions from the isovector Ml 
matrix element (MIV) between the initial ^Sq (T = 1) and the final deuteron (T = 0) 
state, along with the isoscalar Ml matrix element (MIS) and the isoscalar E2 (E2S) matrix 
element between the initial ^Si (T = 0) and the final deuteron '^Si—'^Di states.!^ Since MIV 
is by far the largest amplitude (see below), the spin-averaged cross section (Tunpoiii^P -^ d'y) 
is totally dominated by MIV. Meanwhile, since the initial ^Sq state has J = 0, the MIV 
cannot yield spin-dependent effects, whereas MIS and E2S can. This means that, as first 
noted by Breit and RustgiH, the spin-dependent observables in (ffl) are sensitive to small 
isoscalar matrix elements. To explain the significance of this feature in the context of 
nuclear effective field theory (EFT), we first give a brief discussion of the "chiral filter 
hypothesis" [|]. 

The chiral filter hypothesis based on current algebra was proposed in 1978 prior 
to the advent of modern nuclear EFT. It dictates that corrections to the single-particle 
transitions in nuclei for the isovector Ml operator and the weak axial charge operator 
should be dominated by one-soft-pion-exchange two-body currents. It was shown in ^ 
that this dominance is naturally explained by chiral perturbation theory; the soft-pion 
exchange diagrams are leading-order two-body diagrams in chiral counting. A thorough 
nuclear EFT calculation of MIV carried out in [^ gave a precise, parameter-free estimate 
of CFunpoh which was in perfect agreement with the experimental value |0]; furthermore, it 
was confirmed that the major corrections to the one-body (single-nucleon) current in MIV 
indeed come from one-soft-pion-exchange. For the axial charge also, its large enhancement 
due to the soft-pion exchange two-body currents is by now well established |^. Thus, for 
the observables "protected" by the chiral filter (i.e., observables essentially determined by 
leading chiral-order terms), nuclear EFT has been highly successful. 

The question now is whether we can profitably apply nuclear EFT to observables "un- 
protected" by the chiral filter. With leading chiral-order contributions vanishing, one must 
consider higher order diagrams, a task which can quickly become formidable. The above- 
mentioned isoscalar matrix elements, MIS and E2S, are examples of the "unprotected" 
observable. As is well known, the one-body contribution of MIS is highly suppressed due 
to the orthogonality between the initial ^Si and the final deuteron state. The soft-pion ex- 

*^ The reason why E2S needs to be considered wiU be explained later in the text. We shall also discuss 
later the El contribution recently considered by Chen, Rupak and Savage jj. 



change is also suppressed, there being no isoscalar B^irNN vertex in the leading order chiral 
Lagrangian, where B^ is the external isoscalar field that couples to the baryonic current. 
Due to this double suppression, the size of MIS becomes even comparable (see below) to 
that of E2S, which is a higher-order multipole and hence, in normal circumstances, can be 
ignored. This situation suggests that we must go up to an unusually high chiral order before 
getting sensible estimates of the isoscalar matrix elements that govern the spin observables 
in (HI). However, we show in this paper (confirming the previous results reported in Q) that 
nuclear EFT allows us to make a systematic and reasonably reliable estimation of MIS and 
E2S. 

By the very nature of EFT, an effective Lagrangian involves, at each order of power 
counting, a certain number of parameters. These parameters are in principle calculable 
from first principles for a given scale {e.g., lattice QCD [^) but, in practice, they have to be 
fixed by experimental data. Once the parameters are so fixed, the Lagrangian can then be 
used to make predictions for other processes governed by the same parameters. However, it 
often happens that after fixing the parameters there is not much room left for prediction. A 
major point of our work is to demonstrate that one can make a genuine EFT prediction even 
for the highly suppressed MIS and E2S. Our bona fide prediction not biased by pre-existing 
numbers can be tested soon by the forthcoming results from the ILL experiment Q. 

2 Polarization observables 

Adopting the convention of Miiller et al [^ , we write the transition amplitude as 

{^dnCkA)\T\^np) = xliM{k,X)xnp (2) 

with 



(3) 



Mik,X) = V47r V;. i{kxeX)-{ai-a2)MlY 

■ n ^*N ,^ ^ ,M1S _ p^ ,, ^ ^^ .,,E2S 
-i(k X e;^) • (ax -F 02) -^=- + (cri • ka2 • Ca + ^^2 • ka\ ■ e^j 



where Vn is the velocity of the projectile neutron. As is the deuteron normalization factor 
As ~ 0.8850 fm~ ' , and Xd iXnp) denotes the spin wave function of the final deuteron 
(initial np) state. The emitted photon is characterized by the unit momentum vector k, 
the energy to and the helicity A, and its polarization vector is denoted by e^ = ix{k). The 
amplitudes, MIV, MIS and E2S, represent the isovector Ml, isoscalar Ml and isoscalar E2 
contributions, respectively.^ These quantities are defined in such a manner that they all 



have the dimension of length, and the cross section for the unpolarized np system takes the 
form 

a^^po, = |M1V| V |M1S| V |E2S|^ (4) 
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These amplitudes are real at threshold. 



As we shall see below, the isoscalar terms (|M1S| and |E2S| ) are strongly suppressed 
relative to |M1V| — approximately by a factor of ~ O(10~^) — so the unpolarized cross 
section is practically unaffected by the isoscalar terms. As mentioned, the isovector Ml 
amplitude was calculated [^ very accurately up to 0{Q^) relative to the single-particle 
operator. The result expressed in terms of MIV is: MIV = 5.78 it 0.03 fm, which should 
be compared to the empirical value \ (y'^^poi ~ ^•'''^l '^ 0.004 fm. In this paper, therefore, 
we will focus on the isoscalar amplitudes. 

The isoscalar matrix elements are given by 



MIS 




E2S 


- J'^^*'-'^^''^*'-'^ 


fl = 


- (fxxx Jem(^), 


Q'^ = 


fd^x{3x'x^ -6'^x'^)J^y^{x), 



(5) 
with 



(6) 

where J^^{x) is the EM current. The spin-triplet initial np and the final deuteron wave- 
functions, iV't^) and iV'rf^)) are given as 



(rli^t)- ^ 



'Anr 
with the boundary conditions 



ut,d{r) H y^wt^ir) 



XuXoo (7) 



lim Udir) = Age '^''■'^ , lim ut{r) = r — at , (8) 

r— >oo r—*oo 

and 

,. Wd{r) wt{r) 

hm — — - = r]d , hm -—- = , (9) 

r^oo Ud(r) r^oo Ut{r) 



where jd = y/^dJ^N with Bd — 2.224575 MeV the binding energy and rriN — 939 MeV the 
nucleon mass; rjd ^ 0.0250 is the asymptotic D/S ratio of the deuteron; at ^ 5.4192 fm is 
the ^5i np scattering length, Suir) = 3ai • f a2 • f — ai ■ a2, and x (C) is the spin (isospin) 
wavefunctions. 

As mentioned, the spin-dependent observables are sensitive to the isoscalar matrix 
elements. Let I\{9) be the angular distribution of photons with helicity A = ±1, where 9 
is the angle between k (direction of photon emission) and a quantization axis of nucleon 
polarization. For polarized neutrons with polarization P„ incident on unpolarized protons, 
the photon circular polarization P^ is defined by 

J+i(0°)-/-i(0°) 
/+i(0°) + /_i(0°)' 



p. - 7ZJ. rZi ■ (10) 



At threshold^, P^ is given by 



J- ^ '^Ml + '^E2 

where we have defined the ratios 

_ MIS _ E2S , , 

^^^ = M1V' ^^^ = M1V- (^^) 

When both protons and neutrons are polarized (along a common quantization axis) with 
polarizations Pp and Pn, respectively, the anisotropy r/^ is defined by 

^ /(90°)-/(0°) 
^7 - /(90o)+/(oo)' ^'"^^ 

where /(0) = /+i(0) + /_i(0) is the angular distribution of total photon intensity (regardless 
of their helicities). At threshold, we have 

„ = pP ^Ml + ^E2 - 6^Ml7^E2 .-^^x 

'"' "^ 4{i-pP) + {4 + pP){nl,^ + nl^) + 2pPnMinE2 ^ ' 

where 

pP = Pp-Pn. (15) 

We note that P^ is linear in TZ^s in the leading order. By contrast, unless pP ~ 1, the 
anisotropy r]^ is quadratic in TZ^s and hence highly suppressed. For pP ~ 1 (nearly perfect 
polarization), tj^ becomes 0(1) in powers of 7^'s and can be of a substantial magnitude. 
This feature suggests the possibility of obtaining two independent constraints on 7^'s, one 
from Py and the other from rjy, which would allow extraction of the individual values of T^mi 
and TZe2- There exist data on P^ due to Bazhenov et al. |l^, and the ILL experiment |l|] 
is expected to give information on r]^. 

Our argument above is based on the strict threshold kinematics. A few remarks are 
in order here. When the velocity of the projectile neutron Vn goes to zero, the cross section 
diverges as v~^ and consequently the matrix elements as ?;n ^- However, the quantity that is 
actually measured in experiment is the yield, the flux times the cross section. Since the flux 
is proportional to w„, VnCy near threshold is constant with good accuracy, and one measures 
the thermal average of this f„-independent yield. It is customary to translate the measured 
yield into the cross sections using a fixed velocity, Vn = 2200 m/s, which corresponds to the 
average neutron velocity at room temperature. We adopt this convention throughout this 
paper, and we mean by "at threshold" that we are neglecting any higher order corrections 
in Vn- Chen et al.[^ pointed out that at this velocity the isovector El matrix element (we 
denote it by ElV) connecting initial ^Pj waves with the deuteron state can compete with 
MIS and E2S. They calculated P^ including the ElV contribution but the corresponding 
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See below for a caveat on the kinematics of the experiments. 



result for rj^ has not been reported. The inclusion of ElV affects the asymmetry A'J^^i^O) 
defined and calculated in H. However, for photons emitted in the 9 = direction, there 
is no contribution from ElV. This means that Pj, eg. ([lO|) , is not influenced at all by the 
ElV contribution. Meanwhile, according to Q, rj^ arising from ElV has a very different 
pP dependence from that arising from MIS and E2S. In principle, therefore, it should be 
possible to distinguish experimentally these two ry^'s of different origins. In the present 
work limited to the threshold kinematics, we do not discuss the contribution of ElV.p^ 



3 ChPT calculation 

The chiral counting rules relevant to electroweak vertices in two-nucleon systems are 
well known by now (see, e.g., [y]). The chiral index u oi a given irreducible diagram is given 
by 

u = l-2C + 2L + J2^u Ui = di + ei + 'j-2, (16) 

i 

where C is the number of disjoint pieces (C = 2 for one-body currents and C = 1 for two- 
body currents), and L the number of loops; dj, ej and rij are, respectively, the numbers of 
derivatives/rriTr's, external fields and nucleon lines at the i-th. vertex. Since chiral symmetry 
guarantees fj > 0, eq.(|l6|) implies that the one-body EM current Jfg starts at 0{Q"^) 
while the two-body current J^g starts at 0{Q~^). As spelled out in Q, there are additional 
suppression factors of either a kinematical or symmetry origin. For instance, the space-part 
of the one-body current is suppressed by one power in Q relative to its time part. The 
detailed counting rules for the two-body isoscalar EM current are rather complicated and 
will be discussed later. 

At present, there are essentially two "alternative" ways of power counting in setting up 



an effective field theory for two-nucleon systems. One is the original Weinberg scheme |11|, 
in which the leading four- Fermi contact interaction and a pion-exchange are treated on the 
same footing in calculating the "irreducible" graphs for a potential that is to be iterated 
to all orders in the "reducible" channel. The power counting is done only for irreducible 
vertices. The other is the "power divergence subtraction" (PDS) scheme ||l2], in which only 
the leading (nonderivative) four-Fermi contact interaction is iterated to all orders, with the 
higher-order contact interactions and the pion exchange treated perturbatively. While the 
PDS scheme is perhaps more systematic in power counting, we believe that the Weinberg 
scheme is not only consistent with the strategy of EFT but also, in the sense developed 
below, more flexible and predictive with possible errors due to potential inconsistency in 
power counting generically suppressed. In our work we adopt the Weinberg scheme. In 
this framework, the power counting for electroweak transition amplitudes simply reduces 



** In connection to the above-mentioned thermal averaging, we note that the ratio EIV/MIV is propor- 
tional to Vn, while 75.M1 and 7?.e2 are «„- independent. 



to chiral counting in the irreducible vertex for the current, as the current appears only 
once in the graphs. We are then allowed to separate the current matrix element into the 
one-body (single-particle) and two-body (exchange-current) terms. This framework has the 
advantage that it can be straightforwardly applied to n-body systems for n > 2. 

Iterating the irreducible vertex to all orders in the reducible channel is equivalent 
to solving the Schrodinger equation with a corresponding potential. This then suggests 
that we use, in calculating nuclear responses to slowly varying electroweak probes, those 
wave functions computed with so-called phenomenological "realistic potentials". A prime 



example of these realistic potentials is the Argonne vis potential |13| (called in short Avl8). 
Thus in this approach, referred to as the hybrid approach, the one-body and two-body 
transition operators derived from EFT are sandwiched between the initial and final two- 
nucleon wavefunctions distorted by a phenomenological N-N potential. This procedure of 
mapping effective field theory to realistic wave functions has recently been justified by means 



of a cutoff regularization |14]. A similar argument in support of such a hybrid procedure 



has been presented by van Kolck |15]. 



Previously we adopted this scheme to study the unpolarized p{n,^)d reaction for the 
thermal neutron. The cross section for this process was computed to next-to-next-to-leading 
order (NNLO) in the above-described chiral counting scheme for the transition vertex Q, 
and with the use of the Avl8 wave functions. The resulting theoretical cross section, 
^wnvoi ~ ^^^ '^ 2 ™b, agrees very well with the experimental value crfl^ , = 334.2 it 0.5 fm 



'unpol - -J-J-* -L o iiiu, aj^icco vci^. wcii wiuii tiic ci^jjciiiiiciiudi vciiuc u ^^^^^ 

10]' '^unpoi consists of the leading 1-body contribution, 305.6 mb, and the 2-body exchange- 
current contribution, which in conformity with the chiral filter is dominated by the soft 
one-pion exchange. 

As shown in [|lj], the single-particle matrix element has negligible uncertainty, and 
therefore the theoretical error is entirely attributable to the uncertainty in the exchange- 
current operator associated with the short-distance part of the interactions. This short- 
range behavior is not fully controlled in the formulation of [^. The short-distance part 
introduces a scale and renormalization-scheme dependence, and only those results that are 
not sensitive to this uncertainty can be trusted. In a framework that uses realistic wave 
functions, the short-distance scale is set by a (momentum) cutoff proportional to r~^, where 
Tc is the "hard-core radius" that removes the part of the wave function for r ^ Tc 7^ 0. The 
net effect of this cutoff is that in addition to cutting the radial integrals in the coordinate 
space, it removes all zero-range terms in the current operator including zero-ranged counter 
terms. After the four-fermion counter term is removed by the hard core, there are no more 
parameters in the theory ^. This procedure, familiar to nuclear physicists, may be viewed 



as a form of exploiting the scheme dependence. We shall refer to this procedure as the 
hard-core cutoff scheme (HCCS). HCCS can be justified for the process in question by 



*^See appendix B in the second reference of H] for the zero-ranged counter term in question. In the 
framework of |ll3] , which avoids scheme dependence at the expense of predictivity, this counter term plays a 
much more pronounced role, rendering a bona fide prediction infeasible [[L6| . 



using a (momentum) cutoff ~ J'cT^ and sli owing that tlie counter-term contribution "killed" 
by the hard core lies within that uncertainty range .p^ The 1% theoretical error assigned to 
'^unpoi s-bove represents uncertainty in this cutoff procedure. We will see below that, when 
the chiral filter does not apply, the removal of the zero-range counter terms by the simple 
hard-core cutoff may be much less reliable. 

3.1 One-Body Contributions 

Either starting from the well-known EM form factors of the nucleon, or carrying out 
an explicit ChPT calculation that reproduces the form factors, we arrive at the following 
expressions for the one-body operators: 



fJ'iB{ri,r2) 



2mN 



E 

j=l,2 






l + ri(a^f^ 



P^^^-Pil+oitr 



4m% 



m 



N , 



Q?B(n,r2) 



^E 

1=1,2 



{3ry, - 6'''rf) i±^ (l + ^' 



+ 



(2/ig-l) + (2/.y-l)rf ^^,, ^ ^ 
4m% 



2 



2m% 



pfi 



m 



N , 



with 



(17) 



m 



ab 



-rtip^ X aO' + ^rtip^ x a,f - 5'^'n ■ (pl x <?,), 



(18) 



where //s ~ 0.87981 and /iy ~ 4.70589. The matrix elements (^) can then be obtained by a 
straightforward calculation. For example, the leading order one-body matrix elements are 
given by 



MlSl^^ 



^/2( 



eu;2 



2Jv^mN Jo 



, 3 ...^, f ( \ ( \ Wd{r)wt{r) 
dr -Wd{r)wt{r) + fis lud{r)ut{r) 



E2Sl'^ 



euj2 



20^fU, 



n JO 



drr 



(19) 



Ud{r)wt{r) + Wd{r)ut{r) _ Wd{r)wt{r) 
2 V8 

As mentioned earlier, the one-body MIS is suppressed by the orthogonality 

dr [ud{r)ut{r) + Wd{r)wt{r)] = 0. (20) 

Using the Avl8 wavefunctions, the numerical results for the sum of the LO and NLO are 
given as 

MlSiB(fm) = (-4.192 - 0.105) x 10"^ = -4.297 x 10-^ 

E2SiB(fm) = (1.401 - 0.007) x lO^^ = 1.394 x 10~^ (21) 

*" Since our calculations are carried out in coordinate space, re is related to the hard-core radius used in 
standard nuclear physics calculations. 



We see that the size of the one-body MIS is about 10^ times smaller than MIV, a suppression 
which can be traced to the above-mentioned orthogonality and the smallness of fis as 
compared to ^y. 

3.2 Two-body contributions 



'N/V^- 




Figure 1: Generic diagrams for the two-body isoscalar current. The solid circles include counter- 
term insertions and (one-particle irreducible) loop corrections. The wiggly line stands for the external 
field (current) and the dashed line the pion. One-loop corrections for the pion propagator and the 
ttTVTV vertex need to be included at the same order. 

Generic diagrams for the two-body isoscalar current are depicted in Fig. 1. All the 
propagators and the vertices therein should be understood as the renormalized quantities 
that include all the loops and counter term insertions to a given order. The vertices with 
non-trivial renormalization are represented by big filled circles. The one- loop graphs relevant 
to the filled circle in Fig. la (Fig. 16) are drawn in Fig. 2 (Fig. 3). 



':n 



fall 



'> 



fa2l 



,<^ 



.;V 



(51) 



V ^ 



(52) 



Figure 2: One-loop graphs that contribute to the B^ttNN vertex, B'^ representing the external 
isoscalar field. These diagrams give rise to corrections of 0{Q'^) and higher orders to the leading- 
order one-body term. 

We now proceed to discuss in detail the counting rules for the two-body contributions. 
For easy reference, we let N"LO stand for contributions of order ©(Q") relative to the 
leading one-body term. For the chiral-filter protected MIV, since the isovector 'jttNN 
vertex with i/j = enters at the blob of Fig. la, the leading correction appears at N^LO 
and the one-loops and the counter-terms appear at N^LO. For the isoscalar current, not 
protected by the chiral filter, the counting rules are quite different. The B^ttNN coupling 



/\/>y\;i 



'>y\/\j 



\ / 

\ / 

X 

/ \ 



(5) 




Figure 3: One-loop graphs contributing to the two-body baryonic currents. They come at Oiff^) 
and higher orders relative to the LO one-body term. All possible insertions on the external line are 
understood. 

is highly suppressed by chiral symmetry, and its first non-vanishing term appears at 

ViiBixNN) > 2, UiiB^nNN) > 3. (22) 

The same suppression factors appear in the contact B'^NNNN vertex. As a result, leading 
corrections for the MIS are N^LO, which arise from a chiral Lagrangian [1^, g] of the form 



Co 



9l 



Aii>al3 



87r2 



m. 



{daBji - dpBa) Bj,v^ iA^B^ 



ig^id^B, - d,B^)B,[S^, S^,]B^B^B^ 



(23) 



where B^ is the nucleon field; the definitions of the four- velocity v^^, the spin-operator S!^ 
and the covariant pion-derivative A^ can be found in Q. Although the second term contains 
an unknown constant (74, the first term is free from an additional parameter, for it is given 
by the Wess-Zumino term (related to the Adler-Bell-Jackiw anomaly) the strength of which 
can be determined from the KSRF relation. There is also a Lagrangian consisting of 1/m 
corrections that we shall refer to as the "fixed-term" Lagrangian in which all the coefficents 
are uniquely determined |17], 



C 



fixed 



B,, 



9A 

TUN 



{S^ ■D,vA} 



9A 



{iS^ ■ AD^ + 2iSy D A-D + h.c.) + 



B, (24) 



where qa ~ 1-26 is the axial coupling constant and the ellipses denote terms irrelevant to 
our present study. The chiral index of the first term is f, = 1, while all the other terms 
have Vi = 2. It is to be noted, however, that the i^i = 1 term is proportional to the energy- 
transfer,^ which is of order of Q'^/mj\^. Consequently all the terms in eq.(^4|) give N^LO 
contributions to the Ml operators. There is no counter-term contribution to E2S up to 
N'^LO with which we are concerned. 

The loop contributions enter at N^LO and are finite. Usually, counter-term contribu- 
tions enter - regardless of whether the loops are divergent or not - at the same order as the 

*^ We have replaced the energy transfer v ■ qt — pf — p^ by ■:^^ ^ , where the subscript i is the 

particle index. 
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loop contributions. In our case, however, there is no counter term contribution. Additional 
parameters characterizing the nuclear potential appear at the four-fermion vertex in Fig. 3c, 
through a Lagrangian of the form 



1 



C-o = -i:Y.CA{BTABy 



(25) 



Here F^'s stand for non-derivative operators, and Ca are constant coefficients. Due to the 
Fermi-Dirac statistics, only two of the four possible Ca are independent. 

Summing all the two-body contributions up to N'^LO and Fourier-transforming the 
resulting expressions into coordinate space, we obtain the following isoscalar Ml and E2 
operators: 



fj'2B{r) 



2miy 



UL,Su{f)}cT{r) + hc2{r) 



+ [3rf.S-S)(f,sCUr) + -C^{r] 



+S{fisC4{r)-^C^{r: 



1 

-( 

6 



-?fixod / 



+M?r(r), 



Q'lirl 



(3rV^' - 6'^r^ 



Si2{r)Cl (r) + C2{r) 



(26) 
(27) 



where S = 2(^1 + ^2) and L is the orbital angular momentum operator for the relative 
motion. We have separated the contribution from the fixed-term Lagrangian (V^), 



/^?r(r) 



e TI-T2 mlg\ 
2mN 16mAr/2 



Im 



^ll^L 



yoim^r) -5^{m.^r) 



+ \3£{fr-p)--£ip)-^ J 



y2{m^r) 



(28) 



with 



nx) 



^ijk^j^tkl^tixr 



(29) 



where X is either ^ or ff-p; yo{x) = ^^ and y2{x) = (l + f + -^jyoix) are Yukawa 

functions, and p= — |(V — V) is the momentum operator with the understanding that the 
derivatives act only on wavefunctions. The C's are given as 



Clir) 
CUr) 



Nim% e" 



3 Airr 
2mN 



1 + ^^ + 



+ 



N2,ml 



rrij^r m'^r^ J 247rr Jq (zz)"^ 



dz 



K^ix) 



Hsm^, 



iVwz - iVi 



m 



47rr 



1 + 



+ 



m.,Tr mlr^ 



Nsml 



dz 



Sirr Jo (zz)'- 



Ksix) 



+- 



ml 



dz 



27rr Jo zz 
N2m^ f^ dz 



-—^K,{x) + -^xK,{x) 
3 o 



vrr 



zz 



K^ix) 



(30) 
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and 

C2{r) 



(Ji • (72 



Nim'l e 



3 ^—m-irt' 



Noml /"i dz 



Airr 



vrr Jo zz 



-Ksix) + Vi6if) 



Nsm^ f^ f'd^ [g^g(^) ^ 8^^(^) ^ xKo(x)] 



C'4(r) 



247rr [ Jo zz 
- (3 + 2ai • ^2) / 

JO 

' 2?Ti/v 

-Nwz - Ni 

N-sml (I f^ dz 



1 dz K3{x) + 7Ki{x) 



(2:2:) 



m. 



+ 



^ Airr Svrr Jo zz 



—K3{x) + 



2?n,Ar 



(A^wz + 94) + ^2 



247rr [2 Jo zz 



[9K3{x) + 15iri(x) + 2xKoix)] - 3 



/^5 

(^;^)^ 4 



5(f) 



(31) 



where x 



(1 — z) and Ky{x) is the irregular modified Bessel functions^ of order 



V. The coefficients A^'s are defined by 

iVi = 

N2 = 

Ns = 

Nwz = 



n •T2 



n •r2 



5l 



1287r/4 ' 
128^/4 ' 



(3 + 2ti • r2) 



5l 



1287r/4 



' 1 * ' 2 ^ n 9 po 9 



■487r2/2m2 
and Fl and V2 are two independent coefficients involving Ca, 



(32) 



^1 



-(Ti • CT2 -^— + ^^^ 2^ (^TiCT^CArArATiCj^ + (1 ^ 2)j , 



5^1 + a2) 1/2 ^ iVim.(<?i + <?2) + ^^ E (^1^1 {CAr^r^, ai} rfaj + (1 ^ 2)) (33) 



In the above equations, the terms with 34 and A'wz as well as the fixed-terms are due to 
N^LO contributions, while all the other terms come from N^LO contributions. Fortunately, 
the parameters V2 can be absorbed into g^.^ 



94^ 94 = 94 + N^z + 7, ^2 . 



(35) 



The Vi term, which in naive counting is N'^LO, is further suppressed and plays little role. 
The reason is not hard to understand: Vi appears only in C2(r'), which is accompanied 



** We follow the definition of the Bessel functions given in |lq] . 

*^ Similarly, A'wz can be renormalized to absorb some N*LO contribution, 



AT _. N' - M ^^"^^ 

JVWZ -^ J Vwz = J vwz — 



2771 J\ 



•iVi 



(34) 



Unlike the 34 case, however, this does not reduce the number of unknown parameters, since the value of TVi 
is already known. 
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either by L (for MIS) or by r^ (for E2S). As a consequence, the Vi term contribution has 
the form 



Vi / d^rr^'S^r), (36) 



where n = 4 (for MIS) and n = 2 (for E2S), and eq.(36) is obviously zero. Thus we are 
left with only one parameter g'^. In the next section we will show that we can fix g'^ by 
taking the magnetic moment of the deuteron as input. This will enable us to make a totally 
parameter-free prediction for the ILL experiment. 

4 Renormalization 

The two-body contribution to the MIS, denoted by MIS2B) can be expressed as 

MIS2B = dr J^t{r) —^ 54 / d-'r — - — ^ <^^^ (^)> (37) 



where J-tir) consists of C(r)'s [see eqs.(pC)|, 31)] and the final and initial radial wavefunctions. 



The first integral diverges since J-t{r) contains a piece which behaves like l/r"^ for r ^ 0. 
If the hard-core cutoff scheme (HCCS) defined earlier can be used, this divergence simply 
disappears, as the hard-core kills any contributions from r < Vc {vc = hard-core radius). 
A hard-core would remove altogether the second integral in eq.(p7|). HCCS was applied 
quite successfully to chiral- filter-protected quantities, e.g., MIV [^. The numerical results 
in [^ exhibit only slight rc-dependence, a feature that assures the predictivity of HCCS. 
For processes not protected by the chiral-filter, however, there is no reason to believe that 
HCCS continues to be satisfactory. This is because two-body corrections for these processes 
come from "short-ranged physics," and therefore imposing a hard-core will throw out the 
main part of the relevant physics. For the process (2) this problem becomes even more 
serious because the strong suppression of the one-body MIS in our case means that even 
small changes in the two-body contributions can have a significant influence on the total 
transition amplitude. In particular, since we are calculating up to N LO, it is not justifiable 
to neglect by fiat the g'^ term, which is N^LO. 

For systematic evaluation of the short-ranged contributions along with long-ranged 
contributions, a transparent method would be to impose a momentum cutoff A ~ r~^ in 
performing Fourier transformation. This scheme would allow us to include in a well-defined 
manner the contributions from r < r^ whereas they are simply thrown away in HCCS. 
Technically, however, it is cumbersome to introduce a cutoff in Fourier transformation. 
Fortunately, for our present purposes, we can use a convenient scheme equivalent to, but 
simpler than, the momentum cutoff method. For this scheme, we first note that, the large-r 
contributions being scheme- independent, only the contributions r < re involve scheme- 
dependence. Secondly, while J-t{r) contains both the D-wave and S'-wave radial functions, 
the terms involving the D-wave have higher powers of r and hence receive only minor 
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contributions from the small r region. Then the only quantity that is sensitive to the 
renormalization scheme is the small-r contribution of the term containing Ud{r) and Ut{r). 
For small r, however, these radial functions can be approximated as ^^(r) ~ u^(0)r and 
ut{r) ~ ti^(0)r. The small-r contribution (0 < r < re) of the first integral of eq.(|37|) then 
can be written as 

[ 'drJ^t^u20H{0)A (38) 

Jo 

where A is an rc-dependent but wavefunction-independent constant. This contribution 
can be absorbed into the redefintion of the parameter g'^; the resulting g'^ will become rc- 
dependent. The net effect of this procedure is to remove the r < Tc contribution while 
retaining the Tc-dependent g'^ term. Operationally, we can achieve the same goal by replac- 
ing the delta function attached to g'^ with the delta-shell function. 

For a given value of Tc we can fix g'^ by fitting the deuteron magnetic moment. We shall 
call this scheme the modified hard-core cutoff scheme (MHCCS). A regularization scheme 
similar to MHCCS has often been used in the literature (see, for example, p!9|]). Hereafter 
we drop the prime on g'^ for simplicity. 

In MHCCS, eq.(|37|) is rewritten as 

M1S2B = rdrFtir) _ ^/^ g MrMr,) ^ ^IS^'^^ + MlSlg". (40) 

We determine 54 through the renormalization condition for the deuteron magnetic moment, 

^J'd,2B = J drj'dir) + 2mNg4 -^^ = Z^?'' " ^rf,iB> (41) 

1 47rr^ 
54 = 7^ 27^ A/Xrf, (42) 

A^,^^r-/^rf,iB-/^d!rB^ (43) 

For re = 0.01, 0.2, 0.4, 0.6 and 0.8 fm, we find A^^ = 0.024, 0.020, 0.022, 0.023 and 0.023 
(in units of the nuclear magneton); the corresponding values of 54 are: 54(fm ) = 5.06, 
2.24, 0.73, 0.31 and 0.20. In Table |l] we give the renormalized MIS and its breakdown 
into the various individual contributions, for different values of re spanning a wide range, 
0.01 fm < re < 0.8 fm. We see that, while the finite part exhibits some rc-dependence, the 
total MIS is completely insensitive to the value of re, owing to the renormalization to fit 
the magnetic moment of the deuteron. Our results can be summarized as 

MIS = (-2.85 ±0.01) X 10"^ fm, (44) 
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Tc (fm) 


0.01 


0.2 


0.4 


0.6 


0.8 


M1S^^(N3L0) 


-3.384 


-3.352 


-3.197 


-2.789 


-2.159 


MlSff^CN^LO) 


-0.938 


-0.941 


-0.920 


-0.823 


-0.643 


MlSfg'*^(N3L0) 


-4.322 


-4.293 


-4.117 


-3.611 


-2.802 


MlS|f*^(N4L0) 


2.658 


3.095 


2.709 


2.049 


1.327 


MlSff*^ 


-1.664 


-1.198 


-1.408 


-1.562 


-1.475 


MlSiB + MlSff *" 


-5.961 


-5.494 


-5.704 


-5.859 


-5.772 


MIS 


-2.849 


-2.850 


-2.852 


-2.856 


-2.861 



Table 1: Individual contributions to MIS in unit of IQ--"^ fm. MISIb'*" is the sum of MlSfg"''(N3L0) 
and M1S2b*°(N4LO), while MlSfg^'^lN^LO) is the sum of the contributions from the WZ term, 

We also list the results in the hard-core 
cutoff scheme (HCCS) (without the 54 contribution), MISib + MlSfe *^ with MlSm = -4.297 x 
10^'^ fm. The last row gives the total MIS resulting from the MHCCS regularization. 



M1S^^(N3L0), and the "fixed-term," MlSfg'='^(N3L0). 



where the error bar stands for the rc-dependence. It is highly noteworthy that the Tc- 
independence holds even when (74 changes by a factor of ~ 26. 

We also observe that the individual contributions coming from N^LO and N^LO are 
of the same order as the LO one-body contribution. This feature can be traced to the afore- 
mentioned facts that the LO one-body contribution is suppressed by the orthogonality of 
the wavefunctions, and that the N^LO is suppressed due to the smallness oi gp. About half 
of the N'^LO contributions are cancelled by the N^LO contributions. In obtaining the above 
numerical results, we used the KSRF value, gp ~ 5.85, for the pNN coupling constant 
featuring in the Wess-Zumino term. If instead we use a more conventional value, gp ~ 5.25, 
then the resonance-exchange contribution decreases, leading to a larger cancellation between 
the loop contribution and the N^LO contributions. Even with the use of a unrealistically 
small value of the gp which would cause a complete cancellation among the 2-body terms, 
we have found MIS = —2.91 x 10^^ fm (for the entire range of 0.01 fm < Tc < 0.8 fm), 
which differs only about 2 % from MIS in Table 1. We note that the MEC contribution 
reduces the one-body contribution by about 30 %. 

The two-body contribution to the E2S turns out to be quite small, E2S2B = (0.00 it 
0.01) X 10~^ fm for the whole range of Vc = (0.01 ~ 0.8) fm. Combining the one-body and 
two-body contributions, we obtain 



E2S = (1.40 ± 0.01) X 10"^ fm. 



(45) 



Correspondingly, the two-body contribution to the deuteron quadrupole moment also is 
quite small, Qd,2B = —0.002 fm . This means that the discrepancy between the experimental 
value (0.286 fm ) and the one-body contribution (0.273 fm ) cannot be explained in terms 
of the two-body effects as calculated here. 
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We now discuss briefly the possibility of another renormalization scheme. EarUer, 
we invoked eq.(p^) to eUminate the contribution of the Vi term. However, if we use a 
non-zero value of Tc, the delta- function gets smeared and therefore the Vi term can have 
a finite contribution even though suppressed by powers of Qvc- One then might wonder if 
there exists a possibility to renormalize the deuteron quadrupole moment Q^ as well as the 
magnetic moment by adjusting Vi. This is indeed possible, and we shall refer to this scheme 
as MHCCS*. Since the extension of MHCCS to MHCCS* is straightforward, we simply 
quote the results. For the range of re = (0.01 ~ 0.8) fm, we find MIS = -2.867(1) x 10^^ fm 
and E2S = 1.348(0) x lO'^ fm. Thus MHCCS* and MHCCS give very similar results. We 
have found that, as far as the observables are concerned, the MHCCS* results exhibit even 
a higher degree of insensitivity to Tc than those of MHCCS. Meanwhile, the values of Vi 

and 54 in MHCCS* vary as Vi = (3.1 x 10^ ~ 24.4) fm^ and g^ = (5.05 0.09) fm^, when 

re changes from 0.01 fm to 0.8 fm. While neither Vi nor 5(4 is an observable, the occurrence 
of these large coefficients^^ as well as the huge variation of Vi with respect to Vc is rather 



offending to the naturalness of the theory. So we will not consider this scheme seriously. 

To illustrate the convergence of our chiral expansion, we also present the results that 
are renormalized only up to N^LO. To this order, there are no loops, and E2S does not 
receive any two-particle contribution. The results for MIS again turn out to be very close 
to what we obtained in MHCCS, exhibiting little rc-dependence: For 0.01 fm < Tc < 0.8 fm. 



= (-2.84 ±0.01) X 10"^ fm. 



MIS = (-2.84 ±0.01) X 10"^ fm. (46) 



5 Discussion 



The only "isospin-violating" effect that features in the above consideration is the 
trivial fact that the two-nucleon system can emit both isovector and isoscalar photons. 
We examine here whether any significant additional isospin-violating effects can arise due 
to higher-order electromagnetic (EM) interactions, or isospin impurity in the two-nucleon 
system. Although in normal circumstances these effects give only minor corrections, they 
might be important in calculating highly suppressed amplitudes such as MIS and E2S. 
It is particularly important to assess those isospin-violating effects in the dominant MIV 
which mock the isoscalar amplitude. First, we argue that isospin violation in the initial- 
and final-state interactions cannot influence the observables of our concern. To produce 
any spin-dependent effects, the initial up state must have J 7^ 0, and hence, at threshold, 
should be a ^Si state. If this state could be mixed with ^5*0 (the only other available state 
at threshold), then the dominant MIV transition connecting ^Sq and the deuteron would 
produce a rather signiflcant correction to the genuine MIS. In the two-nucleon systems, 
however, the conservation of J dictates that, even in the presence of isospin-violating inter- 



^^^ Prom eqs.(p3, p3), we can infer that the natural scales of the Vi and 34 are Nim^ ~ 0.09 fm'' and 
TVwz — 0.03 fm'*, respectively. 
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actions, there be no mixing between the ^5*0 and '^Si channels. Next, we need to examine 
radiative corrections to the Ml operator itself. Instead of considering all possible diagrams, 
we study here a typical diagram just to get a rough estimate. One of the diagrams that 
give non-vanishing corrections is obtained from Fig. 36 by replacing the one-pion line in it 
with a photon line. Let ^em be the ratio of the contribution of this radiative correction 
diagram to that of Fig. 36. The EM correction becomes more (less) important for lower 
(higher) momentum transfers. When the momentum transfer is of order of m,r) we find 



„2^2 



■^EM ~ a/ovr — 0.09, where aj^ = tf-tt — 0.08 (/^ ~ 93 MeV). Applied to a specific case 
of MIS, however, this ratio is enhanced by i^v/l^s — 5.3, as the isovector jNN vertex can 
contribute to the isoscalar current due to the isospin breaking EM interaction, and as a 
result the ratio of the EM correction becomes ~ 1/2. Thus the EM correction might be 
substantial compared to the N^LO. We remark, however, that, in MHCCS (and also in 
MHCCS*), this amount of correction would cause negligible changes in the final results. 

To conclude, we recapitulate our predictions in the MHCCS scheme for the ratios of 
the relevant matrix elements and the polarization observables: 

7^Ml = (-0.49 ± 0.01) X 10"^ 

7^E2 = (0.24 ± 0.01) X 10~^ (47) 

and strictly at threshold 

P^ = (-1.04 ± 0.01) X 10"^ (48) 

r]^ = 0.80 ±0.02, (49) 

where the errors represent the variance among MHCCS, MHCCS* and the N^LO MHCCS. 
The Tc-dependences are negligible. The above values of P-y and rj^ are for complete polar- 
ization; that is, |P„,| = 1 for P^ and pP = 1 for r]^. As noted, rj^ decreases rapidly when pP 
becomes smaller; for instance, r]^{pP = 0.96) = 0.62x10"^ and r/^(pP = 0.25) = 0.86x10"^. 
The results with HCCS are rather different ^^, 



HCCS : 7^Ml = (-0.98 ± 0.03) x 10"^ (50) 

P-y = (-1.73 ± 0.05) X 10"^ (51) 

r/^ = 0.53 ± 0.01 , (52) 

where, this time, the errors reflect the rc-dependence. The E2S remains the same as in 
eq.(|3). 

The presently available experimental data for P^ [^] is: 

pexp ^ (_^_5 ^ Q_3) ^ ^Q-3 _ (53) 



*^^These differ from the results given in |8| due to the term M1S2b'^'^(N'^LO) overlooked in the preliminary 
calculation reported there. 
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Although the prediction in the HCCS, (|5l|), gives a P^ closer to (|53|), we regard our MHCCS 
result (47) to be more reliable. The results of Chen, Rupak and Savage 0] for the ratios 



TZmi and 7^e2 computed in the PDS counting scheme of |]T2[ are quite close to ours (^7|). 
Even with the rather generous error estimates in [^ taken into account, this agreement is 
remarkable. ^^ 



While we believe that MHCCS is more reliable than HCCS, we should note that the 
MHCCS and HCCS results for T^mi differ by a factor of ~ 2. To be conservative, we should 
take this difference as uncertainty in the treatment of the short-distance physics involved. 
(7^E2 does not distinguish the two hard-core cutoff schemes.) This sort of uncertainty 
was implied in the chiral filter hypothesis of ref.Q. It has been mentioned in various places 
that, when protected by the chiral filter, the HCCS procedure yields quantitatively accurate 
results. But the other side of the coin of the chiral filter suggests that the HCCS procedure 
is likely to fail when unprotected by the chiral filter, as in the case of the isoscalar matrix 
elements at hand. Our results here indicate how to understand the possible demise of the 
much used hard-core cutoff procedure, a procedure which has long been a standard method 
in nuclear physics, in terms of an effective field theory regularization. What comes out as a 
surprise from this work is that, despite the uncertainty due to the short-distance behavior, 
one can make a reasonably quantitative calculation of the observables that are not protected 
by the chiral filter mechanism. We await with great eagerness the forthcoming results of 
the ILL experiment, which will bring crucial information on this issue. 
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